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Abstract. We consider the hydrogen molecular ion H2 in the presence of a strong 
homogeneous magnetic field. In this regime, the effective Hamiltonian is almost one 
dimensional with a potential energy which looks like a sum of two Dirac delta functions. 
This model is solvable, but not close enough to our exact Hamiltonian for relevant 
strenght of the magnnetic field. However we show that the correct values of the 
equilibrium distance as well as the binding energy of the ground state of the ion, can 
be obtained when incorporating perturbative corrections up to second order. Finally, 
we show that He^'^ exists for sufficiently large magnetic fields. 
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1. Introduction 

Since the discovery of strong magnetic fields in astronomical objects, e.g., in the surface 
of neutron stars, the behaviour of atoms and molecules in these media has become 
a subject of wide interest. Also, high-magnetic-field conditions can be mimicked in 
some semiconductors where a small effective electron mass m* and a large dielectric 
constant e reduce the Coulomb force relative to the magnetic force [11]. For a review 
and references see [21, 23, 11, 22]. In particular, the hydrogen molecular ion H2, has 
been studied under the influence of a strong magnetic field in the last 25 years, see e.g., 
[2, 10, 12, 18, 16, 17, 8, 7, 19, 20, 25, 27, 26] for some references, using variational and 
numerical techniques. This system, as the simplest molecule in nature, can exhibit a 
qualitative behaviour of the internuclear separation as a function of the magnetic field 
that could have applications to the study of more complex molecular systems. 

We consider H^, in a constant magnetic field. We use the Born-Oppenheimer 
approximation, with — e and nie the electric charge and the mass of the electron, 
respectively. In Gaussian units the Hamiltonian is given by 



where 




j^^^ + t:^^^ (2) 



is the Coulomb potential due to the nuclei, and R the internuclear separation. Here 
-f^spin = ■ B and A is the magnetic vector potential. 

In this static situation it is convenient to use the Coulomb gauge V ■ A = 0, and 
the choice A = —\B x r. The magnetic field B is taken parallel to the unit vector z of 
the nuclear axis, which is the energetically most favorable situation [26, 12]. 

The criterium to define the strong magnetic field regime is to ask for a gap between 
Landau levels bigger than an energy of the order of the ionization energy of the hydrogen 
atom. This means 

hwo^ — ^ 27.2 eV = 1 Hartree, (3) 

do 

where ujq = eBQ/m^c = 4.13 x 10^^ seg^^ is the cyclotron frequency and Oq = 
It? /mee^ = 0.53 A, the Bohr radius. That gives us the threshold value of 5, 

^ mgc ^ 2.35 X 10^ G. (4) 
The Hamiltonian can be written as 

H = Hluq 



4f UB\L.UB^\^_^^^iBY 



fp 2 2 \BoJ h 8 \BoJ a^o \BoJ h 

where = xpy — ypx- Choosing atomic units, i.e., nie — h — e — 1, the Hamiltonian 
reads 

= ^ + i/osc - ^L, + BS,-V + ^, (6) 
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where 

H^..-^^^ + ^ix' + y^ (7) 

In this system of units, the magnetic field is measured in units of -Bq = c ~ 137, 
the energy in units of hujQ = 1 Hartree = 27.2 [eV] and the separation R in units of 
ao = 0.53 [A]. Moreover, using cyhndrical coordinates (i.e., x = p cos 93, y = psiiKp and 
z = z), 

= 7f rL^ + u ^fii^ 2Ni - 

2. Spectral decomposition 

In order to diagonalize H, we use that [if, L^] = [H, S^] = [L^, S^] = and the spectral 
decompositions 

= mn^"^) and 5^ = s.n^^^^ (9) 

meX s^=±l/2 

to get a diagonal matrix ii = if (™''*^)n("*) (g) II^*^) where 

H(m,s.) ^ I ^ _^^ + Bs,-V+^, (10) 

and iio™'' = -f 0^o(l^l + 2n + l)ni'"'*. Here Iln™'' denotes the projector of the n*'^ 
eigenvalue of Ho^\ 

The ground state belongs to = — | in the spin sector and to m = in the 
sector, assuming that the result of [1] apphes also to H^. So, we take the term 

it 

where = ^ + if^g^ — ^ — V. Using the projectors of the spectral decomposition 

of iiisc we can write in a matrix form as follows 



(12) 



where Hefj = IIq^'' and 11^ = 1 — Ilcfr, and IIq'^'' the projector over the lowest Landau 
level. The "effective" part of our Hamiltonian, which is the physically relevant, 
hes = neff/i*^°'~2)nefr, Can be written as 

/leff = f - Kff (13) 



00 9 



with 

Kff(^) = B I e-'^Viz,p)pdp, (14) 
Jo 

since lleff projects onto the wave function 

= Vl^-"^. (15) 
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Making the change of variable, u — ^ we get 

POO 

Ves{z) = / e-''V{z,u)du (16) 
Jo 

with 

V{z, u) = i -r + i -r. (17) 

((^-ff + f)' ((-+ff + f)' 

In the sequel we present some necessary mathematical results which will appear in 
full elsewhere [3]. Let C/l be the unitary implementation of the scaling z ^ z/L and 
define the scaled matrix as hi — -^ULh^^^^^U^^. Its singular part reads 

^eff = f - ^V^. (18) 

where Vl — V£ + and 



VS(z) = ; -rdu. (19) 

If we consider the resolvents tl — {hi — ^'^'^ "^s — {hs — ® nefr ® where 



we can show that [3] 



|rL-riff||<-| and \\r^^-rs\\<^ (21) 



with = {h^Q — the resolvent of /ifg, and Ci and C2 are some positive 

constants. The appropriate scaling law for L in terms of the magnetic field is given by 
L = L(B) = 2W{^), where W{x) denotes the Lambert function (for more information 
about this function see [5]). This relationship between L and the magnetic field B is 
the correct scaling law to insure the convergence of the resolvents of the two different 
models which follows from the proof of the key theorem in [3] . Such a resolvent estimate 
implies, in particular, that the ground state of /il, our scaled Hamiltonian, is asymptotic 
to Co, the ground state of hs, when the magnetic field tends to infinity. 

3. Perturbation theory 

In what follows, we will compute the ground state of h^Q perturbatively from the 
asymptotic model hs- 

The ground state energy, cq, of hg is given by 

"0.1 , W(2ae-^^) 
eo^-y with q;o = 1 + ^2^ (22) 

where a = RL/2. Here ao is the solution of the equation aQ = 1 + e~^""°, which can be 
expressed as (22), in terms of the Lambert function W. 
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Since the Hamiltonians hi and hs are invariant under parity and since the ground 
state is even, it is convenient to work from here on in L'^(R^). This being the case, the 
ground state of hs in L^(R+) reads. 



cosh(Q;o-2) z < a 



(23) 



with Ai = dQe^"'^^ A2I2 and A2 = 2/^2(2a + e^-^^o). 
We set 



6{z-a)--VL. 



(24) 



AV = n^ft - ris 

To second order in perturbation theory, the ground state energy of /igg is given by [9] 
62 = eo + tr (Po Ay Po) - tr {P^AVuAVP^) (25) 



where Pq is the projector over -00 and is the corresponding reduced resolvent. 
The free kernel is given by 

1 



-2^\x-y\ 



+ e" 



(26) 



and we use the notation d^Go{x, y; ^) = dGo{x, y; Let Go{x, y) be the free kernel 

evaluated at ^ = cq, and correspondingly we let d^Go{x, y) be d^Go{x, y; ^) evaluated at 
^ = Bq. Using this notation, the kernel of the reduced resolvent at ^ = eo = — c^o/^ is 
given by 

/ Go(a,y) \ 
y d^Go{a,y) J 



G{x, y) = Gq{x, y)+{Go{x, a) d^Go{x, a)) 



A3 Ai 




(27) 



with 



1 dlG^{a,a) 



-1 



(28) 



2 (9^6-0(0, a))'' d^G^{a,ay 

Therefore, we can compute for each B and R the energy of the molecule H2 in a 
magnetic field, now including the scaling and the repulsion energy, using the expression 

1 



E2{B,R)=L\B)e2{B,R) + 



R 



(29) 



The minimum of E2{B, R) as a function of R, determines the equilibrium distance 
between the nuclei as well as the binding energy, which is defined as this minimum. We 
set E2 = E2{B, Reg), where Rgg denotes the equilibrium distance. 
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4. Numerical results 

We have computed, with Mathematica [15], the equihbrium distance and binding energy 
of the H2 molecule using second order perturbation theory, for a wide range of magnetic 
fields. Using the high level language of Mathematica, our program to compute E2{R) is 
simply a transcription in this language of the mathematical formula (25), (29) and (31). 
Simple and double integrals are made with the numerical integration of Mathematica 
using the Gauss-Kronrod method option. Plotting R E2{R) we are able to locate 
Reg and finally to compute E2{Req)- In tables I, II, III and IV, we compare our results 
for the binding energy and internuclear separation, respectively, with those found by de 
Melo et al. [16] (using variational techniques), Le Guillou et al. [12], Lai et al. [10] and 
Heyl et al. [13], respectively. In figures 1 and 2 we plot the binding energy, and in figure 
3 we plot the internuclear equilibrium distance, both against the natural scaling of these 
quantities, i.e., for the binding energy and L for the internuclear distance. Finally, 
in figure 4 we see that the product ReqL is monotonically decreasing as B increases. 

5. Stability of homonuclear diatomic molecules with one electron 

So far we have discussed the behavior of the ground state of H2 for large magnetic fields. 
In this section we will consider homonuclear diatomic molecules with one electron but 
with nuclear charge not necessarily one. We are interested on the stability of these 
molecules as a function of Z and B. In particular, we are interested in determining the 
highest value of Z for which a diatomic molecule of this sort may exist, as a function of 
the strength of the magnetic field. We have not try to go beyond realistic values of B, 
i.e. B > 10^^. We will use the same model as the one discussed above but this time with 
nuclear charge Z. In order to take into account the explicit dependence of the model 
on Z we must change a to a{Z) = ZRL/2 and V^{z) to 



and, correspondingly, the second order perturbation theory expression for the energy of 
the system is given by 



Here 62 is given, as before, by (25), but with the appropriate changes in AV. 

In the case of zero magnetic field, it was shown numerically by Hogreve [14], 
that a homonuclear diatomic molecule with one electron exists as long as Z is below 
Zcrit = 1.2367. From our calculations here (sec Table V below) we see that the magnetic 
field enhances the binding properties of this molecular ion. The maximum Z for which 
a homonuclear diatomic molecule with one electron exists, is a monotonically increasing 
function of B. Concerning our notation in Table V, below Z^^ we have a bound state 




(30) 



E2{B,R,Z)^Z^ L\B)e2{B,R,Z) + ^ 



(31) 
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of the molecule. Between Z^^ and Z*^^ we have only a local minimum, i.e., a resonance. 
Above Z'^^ , we have no minimum for the potential energy curve R — > E2{B, R), see (29), 
at finite values of R. 

6. Concluding remarks 

Applying perturbation theory to a solvable model, we have computed the binding energy 
and the internuclear separation for the H2 molecule in the presence of a strong magnetic 
field. Our results are in good agreement with those found in the literature for a wide 
range of the magnetic field. We also give the critical nuclear charge for the stability of 
diatomic molecules with one electron in a strong magnetic field. We remark that for 
magnetic fields about 10^^ gauss , we can find a diatomic molecule with one electron 
with a nuclear charge Z = 2, which is a new atomic system. Finally, we note that one 
advantage of perturbation theory on variational methods is to give the possibility of 
computing numerically a window which contains the exact value of the ground state 
energy E{R) of the effective Hamiltonian hes- In other words to produce lower bounds 
on E{R) whereas variational methods are not good at that. Even using (21) we may 
compute these windows for the ground state energy of the complete Hamiltonian H. 
We hope to come back to this question in a further work. 
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Table 1. Comparative table of data with those ones 
calculated by de Melo et al.[16] (with a 

B (Gauss) Reg (a.u.) -E2 (Hartree) R* (a.u) -E* (Hartree) 



1 X 10^° 


1.494 


1.49 


1.232 


1.42 


5 X IQi" 


0.813 


2.92 


0.736 


2.77 


1 X 10" 


0.632 


3.84 


0.604 


3.64 


5 X IQii 


0.364 


6.97 


0.351 


6.59 


1 X 10^2 


0.291 


8.86 


0.285 


8.38 


5 X IO12 


0.182 


14.90 


0.179 


14.08 


1 X 10" 


0.148 


18.35 


0.149 


17.32 


5 X IQi^ 


0.099 


28.76 


0.104 


26.90 


1 X 10" 


0.084 


34.40 


0.085 


31.86 



Table 2. Comparative table of data with those ones 
calculated by Guillou et al.[12] (with a 

B (Gauss) Req (a.u.) —E2 (Hartree) R*^ (a.u) —E* (Hartree) 



1.175 X lO^" 


1.403 


1.59 


1.358 


1.54 


2.35 X 10^° 


1.073 


2.14 


1.038 


2.06 


3.525 X 10^" 


0.923 


2.53 


0.893 


2.43 


4.7 X 10^" 


0.830 


2.84 


0.803 


2.73 


5.875 X 10^'' 


0.766 


3.11 


0.740 


2.98 


1.175 X 10^1 


0.596 


4.08 


0.578 


3.91 


2.35 X 10" 


0.467 


5.30 


0.455 


5.08 


4.7 X 10" 


0.371 


6.82 


0.362 


6.54 


7.05 X 10" 


0.325 


7.86 


0.318 


7.55 


1.175 X 10^2 


0.276 


9.36 


0.271 


9.01 


2.35 X 10^2 


0.224 


11.75 


0.221 


11.35 


4.7 X 10^2 


0.183 


14.62 


0.181 


14.17 



Table 3. Comparative table of data with those ones 
calculated by Lai et al.[10] (with a -k). 

B (Gauss) Req (a.u.) -E2 (Hartree) R*^ (a.u) -E* (Hartree) 



1 X 10" 


0.632 


3.84 


0.61 


3.67 


5 X 10" 


0.364 


6.97 


0.35 


6.69 


1 X 10^2 


0.291 


8.86 


0.280 


8.53 


2 X 10^2 


0.235 


11.16 


0.230 


10.78 


5 X 10^2 


0.180 


14.90 


0.180 


14.46 


8 X 10^2 


0.158 


17.18 


0.15 


16.71 


1 X 10^3 


0.148 


18.35 


0.15 


17.88 


1 X 10^^ 


0.084 


34.40 


0.085 


33.83 


5 X lO^* 


0.060 


50.60 


0.060 


50.07 
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Table 4. Comparative table of data with those ones 



calculated by Heyl et al.[13] (with a 


B (Gauss) 


Reg (a.u.) 


-E2 (Hartree) i?*^ (a.u) 


-E* (Hartree) 


9.4 X 


0.151 


18.02 


17.52 


2.35 X 1013 


0.119 


23.43 


22.89 


4.7 X 1013 


0.100 


28.29 


27.69 


9.4 X 10^3 


0.085 


33.87 


33.28 


2.35 X 10^4 


0.070 


42.44 


41.73 


4.7 X lO^"* 


0.061 


49.89 


49.14 



Table 5. Data of the stability study of the system 

Z=l Z=1.2 Z=1.4 

B (Gauss) Reg (a.u.) —E2 (Hartree) Rgg (a.u.) —E2 (Hartree) Rgg (a.u.) —E2 (Hartree) 



1 X 10^° 


1.494 


1.49 


1.718 


1.71 






1 X 10" 


0.632 


3.84 


0.712 


4.51 


0.810 


5.19 


1 X 10^2 


0.291 


8.86 


0.318 


10.58 


0.353 


12.27 


1 X 10^3 


0.148 


18.36 


0.157 


22.23 


0168 


25.98 


1 X 10" 


0.084 


34.40 


0.087 


42.33 


0.090 


50.05 



Z=1.6 Z=1.8 Z=2 



B (Gauss) 


Reg (a.u.) 


—E2 (Hartree) 


Reg (a.u.) 


-E2 (Hartree) 


Reg (a.u.) 


—E2 (Hartree) 


1 X IO12 


0.400 


13.99 


0.445 


15.78 






1 X 10" 


0.183 


29.67 


0.202 


33.39 


0.227 


37.26 


1 X 10^^ 


0.095 


57.55 


0.101 


64.88 


0.108 


72.14 



Z=2.2 Z=2.4 


B (Gauss) 


Reg (a.u.) —E2 (Hartree) Reg (a.u.) 


—E2 (Hartree) 


1 X 10^^ 


0.118 79.44 0.131 


86.98 



Table 6. Upper bounds to the critical nuclei charge of stability of the molecule given 



by our model. 



B (Gauss) 






1 X 10^° 


< 1.32 


< 1.58 


1 X 10" 


< 1.55 


< 2.05 


1 X 10^2 


< 1.80 


< 2.60 


1 X 10^3 


< 2.10 


< 3.10 


1 X 10^^ 


< 2.43 


< 3.50 



H2 in a strong magnetic field 



11 



-Binding Energy 
50 - 



X 


Heyl 







Lai 




A 


Melo et 


al. 


□ 


Guiiiou 


et ai 


- 


- Present 





50 60 



Figure 1. Comparison of the binding energies for all magnetic field values. The 
binding energy is plotted in atomics units. The adimensional scaling parameter L is 
calculated as L = 2W(|y^^) and therefore is independent of system of units chosen 
for the magnetic field. 




Figure 2. Detail of the figure 1 
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Figure 3. Comparison of th.e equilibrium internuclear distance, which is plotted in 
atomic units. 
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Figure 4. Behaviour of the product of Re, L 



